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Abstract 

Mechanisms  by  which  trace  amounts  of  CO  and  C02  impurities  in  fuel  may  affect  the  performance  of  direct  hydrogen  polymer-electrolyte  fuel 
cell  stacks  have  been  investigated.  It  is  found  that  the  available  data  on  CO-related  polarization  losses  for  Pt  electrodes  could  be  explained  on  the 
basis  of  CO  adsorption  on  bridge  sites,  if  the  CO  concentration  is  less  than  about  100  ppm,  together  with  electrochemical  oxidation  of  adsorbed 
CO  at  high  overpotentials.  The  literature  data  on  voltage  degradation  due  to  C02  is  consistent  with  CO  production  by  the  reverse  water-gas  shift 
reaction  between  the  gas  phase  C02  and  the  H2  adsorbed  on  active  Pt  sites.  The  effect  of  oxygen  crossover  and  air  bleed  in  “cleaning”  of  poisoned 
sites  could  be  modeled  by  considering  competitive  oxidation  of  adsorbed  CO  and  H  by  gas  phase  02.  A  model  has  been  developed  to  determine  the 
buildup  of  CO  and  C02  impurities  due  to  anode  gas  recycle.  It  indicates  that  depending  on  H2  utilization,  oxygen  crossover  and  current  density, 
anode  gas  recycle  can  enrich  the  recirculating  gas  with  CO  impurity  but  recycle  always  leads  to  buildup  of  C02  in  the  anode  channels.  The 
buildup  of  CO  and  C02  impurities  can  be  controlled  by  purging  a  fraction  of  the  spent  anode  gas.  There  is  an  optimum  purge  fraction  at  which 
the  degradation  in  the  stack  efficiency  is  the  smallest.  At  a  purge  rate  higher  than  the  optimum,  the  stack  efficiency  is  reduced  due  to  excessive 
loss  of  H2  in  purge  gas.  At  a  purge  rate  lower  than  the  optimum,  the  stack  efficiency  is  reduced  due  to  the  decrease  in  cell  voltage  caused  by  the 
excessive  buildup  of  CO  and  C02.  It  is  shown  that  the  poisoning  model  can  be  used  to  determine  the  limits  of  CO  and  C02  impurities  in  fuel  H2  for 
a  specified  maximum  acceptable  degradation  in  cell  voltage  and  stack  efficiency.  The  impurity  limits  are  functions  of  operating  conditions,  such 
as  pressure  and  temperature,  and  stack  design  parameters,  such  as  catalyst  loading  and  membrane  thickness. 

©  2008  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Whereas  Pt  electrodes  in  low-temperature  polymer- 
electrolyte  fuel  cells  (PEFC)  show  very  small  polarization  losses 
for  the  hydrogen  oxidation  reaction,  the  losses  can  be  signifi¬ 
cant  if  the  fuel  H2  is  not  pure  but  contains  small  amounts  of 
CO  [1].  Furthermore,  because  automotive  PEFC  stacks  oper¬ 
ate  at  substantially  less  than  100%  H2  utilization  per  pass,  the 
spent  anode  gas  containing  unconverted  H2  must  be  recycled  to 
the  anode  inlet  and,  thus,  may  enrich  in  CO  and  other  impuri¬ 
ties,  if  they  are  present.  Natural  factors  that  limit  the  buildup 
of  CO  in  anode  channels  include  electrochemical  oxidation  of 
CO  to  CO2  at  high  overpotentials  and  its  reaction  with  O2  that 
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crosses  over  from  the  cathode  channels.  Also,  the  buildup  of 
CO  and  other  impurities  can  be  actively  controlled  by  purg¬ 
ing  a  fraction  of  the  spent  anode  gas  recirculated  to  the  anode 
inlet,  albeit  at  the  expense  of  H2  that  is  lost  with  the  purge 
gas. 

Hydrogen  produced  by  reforming  hydrocarbons  may  contain 
CO2  in  addition  to  CO  and  other  reactive  and  inert  impurities. 
Carbon  dioxide  is  not  adsorbed  as  strongly  on  Pt  as  CO  but 
is  not  completely  inert  either.  The  mechanism  of  CO2  poison¬ 
ing  of  the  Pt  catalyst  may  involve  CO2  reduction  to  CO  by  H2 
either  electrochemically  or  through  the  reverse  water-gas  shift 
reaction. 

The  purpose  of  this  paper  is  to  analyze  the  effects  of  trace 
CO  and  CO2  impurities  in  fuel  H2  on  the  performance  of  the 
PEFC  stacks  for  automotive  applications.  More  specifically,  the 
purpose  is  to  address  the  following  questions  related  to  these 
impurities: 
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Nomenclature 

A  membrane  area 

E  potential  (or  activation  energy) 

/p  purge  fraction 

F  Faraday  constant 

AH  enthalpy  of  desorption 

J  current  density 

k  rate  constant 

M  active  sites 

N  molar  flow  rate 

Nc  catalyst  site  density 

P  pressure 

R  gas  constant  (or  recycle  ratio) 

T  temperature 

V  cell  voltage 

z  charge  transfer 

Greek  letters 

a  anodic  transfer  coefficient 

am  coefficient  of  net  water  transport  across  mem¬ 
brane 

symmetry  factor 
8  thickness 

y)  overpotential 

k  water  content  of  membrane 

0  fractional  coverage 

ifr  permeance 

Subscripts/superscripts 
a  anode  (or  adsorption) 

ae  anode  electrode 

c  cathode 

ce  cathode  electrode 

d  desorption 

e  electronic  (or  electrochemical) 

g  gas 

i  ionic 

m  membrane 

N  Nemst 

p  bipolar  plate 

r  reaction  (or  reverse  reaction) 

ref  reference 

s  activation  (or  surface  reaction) 


(a)  What  are  the  mechanisms  by  which  CO  and  CO2  impurities 
in  fuel  H2  affect  the  performance  of  fuel  cells? 

(b)  What  is  the  effect  of  anode  gas  recycle  on  buildup  of  these 
impurities? 

(c)  How  does  the  buildup  of  CO  and  CO2  depend  on  purge  rate? 

(d)  What  is  the  effect  of  the  buildup  of  impurities  on  cell  volt¬ 
age? 

(e)  What  are  the  impacts  of  purge  and  impurity  buildup  on  stack 
efficiency? 


(f)  What  might  be  the  “acceptable”  levels  of  these  impurities 
in  fuel  H2? 

The  focus  of  this  work  is  on  the  types  of  Pt  electrodes  and 
structures  that  are  being  considered  for  automotive  direct-H2 
fuel  cells.  As  such,  we  have  not  considered  alloys  in  which  a 
second  element,  such  as  Ru,  Sn,  or  Mo,  is  added  to  improve 
the  CO  tolerance  of  the  Pt  catalyst.  These  types  of  alloys  are 
more  suitable  for  fuel  cells  operating  with  reformate  streams 
that  may  contain  impurities  at  higher  concentrations  than  in  fuel¬ 
cell  quality  hydrogen.  In  validating  our  model,  we  have  used 
published  experimental  data  obtained  with  finely  dispersed  Pt 
on  high-surface  area  carbon  supports.  With  some  modifications, 
our  model  may  be  applicable  to  other  electrode  structures,  such 
as  nanostructured  thin  film  catalysts,  and  alloys  in  which  the 
second  (or  third)  element  is  added  to  improve  the  mass  activity 
and  durability  of  Pt. 

2.  Model  formulation 

As  briefly  outlined  in  Refs.  [2,3],  we  solve  a  set  of  equations 
for 

(a)  the  ionic  and  electronic  potential  distributions, 

(b)  current  generation  in  the  catalyst  layers, 

(c)  ionic  and  electronic  current  distributions, 

(d)  water  transport  across  the  polymer-electrolyte  membrane, 

(e)  H2,  O2,  N2  and  H2O  concentrations  across  the  porous  cat¬ 
alyst  and  gas  diffusion  media  layers, 

(f)  species  concentrations  in  the  gas  channels,  and 

(g)  capillary  transport  of  liquid  water  across  the  porous  media. 

From  these  solutions,  we  determine  the  Nernst  potential,  E n, 
and  the  overpotentials,  rj,  in  terms  of  which  the  cell  voltage  V 
can  be  written  as 

V  =  EN  -  rj?  -  r]-f  -rt-nt-vl-nl-  mm  -  rg 

+  rf*  -  rjf  -  rg  (1) 

where  the  subscripts  e,  i,  g,  and  s  denote  the  electronic,  ionic, 
gas-phase,  and  activation  components  of  the  overpotentials,  and 
the  superscripts  ae,  a,  m,  c,  ce,  and  p  represent  the  anode  elec¬ 
trode  (gas  diffusion  layer),  anode  catalyst,  membrane,  cathode 
catalyst,  cathode  electrode  (gas  diffusion  layer),  and  the  bipolar 
plate,  respectively. 

The  method  of  calculating  the  various  overpotentials  in  Eq. 
(1)  has  been  described  in  Refs.  [2,3]  and  will  not  be  repeated. 
Here,  we  discuss  the  development  of  a  method  to  determine  the 
effect  of  CO  and  CO2  impurities  on  the  anode  overpotential,  r}\. 

2.1.  Poisoning  model 

In  our  model,  the  kinetics  of  the  hydrogen  oxidation  reac¬ 
tion  (HOR)  are  described  by  considering  dissociative  adsorption 
of  H2  on  Pt  sites  followed  by  electrochemical  oxidation  of  the 
adsorbed  H2  species. 
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2k):o0co  exp 


(i  -  pyeco 
RT 


H2  +  2M^  2M  -  H 

(R-l) 

kd 

k& 

M  -  H^>M  +  H+  +  e“ 

(R-2) 

Camara  et  al.  [4]  have  shown  that  both  linearly  bonded  and 
bridge-bonded  species  must  be  considered  in  order  to  explain 
the  experimental  data  on  CO  poisoning  of  Pt  electrocatalyst. 
Our  analysis  suggests  that  the  low  CO  concentration  data  can 
be  adequately  explained  by  considering  only  the  bridge-bonded 
sites.  Since  our  interest  is  in  Pt  poisoning  at  low  CO  concen¬ 
trations,  the  model  can  be  simplified  by  considering  associative 
adsorption  of  CO  on  bridge  sites  and  electrochemical  oxidation 
of  bridge-bonded  CO. 

Kco 

CO  +  2M  <=*  M2  —  CO  (R-3) 

Kco 

ke 

M2  -  CO  +  H20-^-2M  +  C02  +  2H+  +  2e~  (R-4) 

Analyses  by  Janssen  [5]  and  others  indicate  that  the  effect  of  CO2 
on  cell  polarization  data  can  be  simulated  by  considering  the 
catalytic  reverse  water-gas  shift  reaction  (RWGS).  Consistent 
with  reactions  (R-3)  and  (R-4),  we  include  the  following  reaction 
between  the  adsorbed  H2  species  and  gaseous  CO2  to  produce 
bridge-bonded  CO. 

krc  o 

C02  +  2M  -  H— >M2  -  CO  +  H20  (R-5) 

For  very  low  CO  concentrations  and  thin  membranes,  enough 
O2  may  cross  over  from  the  cathode  air  to  affect  the  site  balance. 
We  model  the  effect  of  this  internal  oxygen  bleed  by  including 
the  competitive  oxidation  of  the  adsorbed  CO  species  and  H2 
species. 


1  ks 

M2  -  CO  +  -  02  2M  +  C02  (R-6) 

2M  —  H  +  i-02— 4>2M  +  H20  (R-7) 

We  define  6m,  and  6qo  as  the  fractions  of  Pt  sites  that  are, 
respectively,  vacant,  occupied  by  H  and  occupied  by  CO.  We 
assume  that  there  is  only  one  type  of  site  and  that  all  sites  can 
be  bridge-bonded  by  CO.  It  is  further  assumed  that  the  reaction 
order  is  the  same  as  the  molecularity  for  all  species  in  reactions 
(R-l)  through  (R-7),  except  O2  in  reactions  (R-6)  and  (R-7)  for 
which  the  reaction  order  is  taken  as  1  rather  than  1/2.  Analysis 
of  experimental  data  showed  a  better  match  with  theory  if  the 
reaction  order  is  taken  as  1.5  for  M— H  in  reaction  (R-2)  and 
CO2  in  reaction  (R-5).  With  these  assumptions,  the  steady-state 
site  coverage  can  be  determined  by  solving  the  following  set  of 
algebraic  equations: 

^upu2^m  -  -  2/^6^' 5  sinh  ^ 

~  ^co2^co2^  5^h  —  2&h$h^o2  =  0  (2) 


2kCOPCO»M  exP  (  #7  °)  “ 

~  2kQO0co  sinh  +  2k'COl(  P^ok)'5 

-2Cco0coPlh=()  (3) 

In  writing  Eqs.  (2)  and  (3)  we  have  further  assumed  that  the  kinet¬ 
ics  of  the  electrochemical  oxidation  reactions  can  be  described 
by  Butler-Volmer  type  equations.  Also,  a  correction  factor  has 
been  applied  to  CO  adsorption  and  desorption  rates  as  per  the 
Temkin  adsorption  isotherm  [6]. 

Knowing  site  coverages  and  reaction  rates,  the  variation  of  the 
molar  flow  rates  of  H2,  N2,  CO,  CO2, 02  and  H2O  along  the  flow 
direction  can  be  determined  from  the  following  conservation 
equations: 


dNU2 

dA 

7h  VfH2  Pu2 

IF 

-  NcksHd^P^2  -  Nckco2Pco2^u 

(4) 

d  1Vn2 

_  M(Pk  -  P$,2) 

(5) 

dA 

5m 

diVo2 

dA 

—  l/2Nc(kco6coPo2  ku0uPo  ) 

5m 

(6) 

d/Vco 

dA 


Nc 


— ^CO^CO^MexP 


+  kco°co  exp 


(1  -  pyeco 
RT 


(7) 


+  ^CO^2-^o2(^CO2)1'50h)  (8) 

diV^o  =  +  A,c(^o^2(/?oJI.5  +  (9) 


In  the  above  equations,  the  current  density  is  related  to  the 
electrochemical  oxidation  rates  as 

yH  =  2FNckeH0lH5  sinh  (10) 

Jco  =  2FNckecodCo  sinh  (11) 

J=Ju  +  Jco  (12) 


2.2.  Method  of  solution 

An  implicit  finite-difference  scheme  was  used  to  solve  the 
governing  equations.  The  fuel  cell  was  divided  into  5-50  axial 
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Fig.  1.  Anode  gas  recycle  with  partial  purge. 


nodes  and  a  marching  algorithm  was  used  to  solve  the  result¬ 
ing  non-linear  algebraic  equations  from  cell  inlet  to  outlet.  The 
solution  was  iterative  because  of  the  recycle  loop  that  couples 
the  anode  inlet  boundary  conditions:  (1)  to  the  gas  compositions 
of  fuel  H2  (0)  and  the  anode  effluent  (2)  as  shown  in  Fig.  1. 

Nf  =  vf  +  (1  -  fp)Nf\  (13) 

where  j  =  H2,  N2,  CO,  CO2  or  O2. 

A  simple  fixed-point  iteration  method  was  used  to  carry  out 
the  iterations.  Usually,  fewer  than  100  iterations  were  needed  to 
achieve  convergence  to  any  desired  accuracy. 

3.  Model  validation 


Extensive  data  is  available  in  the  literature  that  could  be  used 
to  validate  the  CO/CO2  poisoning  model  presented  in  Section 
2.  We  specifically  selected  data  to  validate  the  models  for  CO 
poisoning  as  functions  of  CO  concentration  and  temperature, 
effect  of  O2  in  relieving  the  CO  poisoning  and  reverse  water  gas 
reaction  between  adsorbed  CO  and  gaseous  O2. 

Table  1  summarizes  our  choice  of  rate  and  model  constants  for 
the  adsorption,  desorption,  electrochemical  oxidation,  RWGS 
and  surface  oxidation  steps  in  reactions  (R-l)  through  (R-7).  The 
constants  in  Table  1  are  given  for  85  °C  temperature  (Tref)  and 
have  been  either  taken  from  literature  or  derived  from  data  fitting. 
The  rate  constants  for  the  adsorption  steps  are  scaled  with  respect 
to  temperature  based  on  the  kinetic  theory  of  gases  applied  to 
molecular  collisions  with  surfaces  and  for  the  desorption  steps 
on  the  basis  of  the  heat  of  desorption. 


0.5 


kt  =  kt,ref  exP 


ah) 

R 


(14) 

(15) 


where  i  denotes  H  in  reaction  (R-l)  and  CO  in  reaction  (R- 
3).  We  have  inferred  a  A//^  of  9.5  kJ mol-1  and  a  A H^0  of 
150  kJ  mol-1  from  the  analysis  of  the  selected  data  which  is  in 
agreement  with  the  data  reported  in  Ref.  [7] . 


Table  1 


Model  constants  derived  from  experimental  data  reported  in  Refs.  [1,4, 8, 9] 


Variable 

Unit 

% 

mol/(mpt  bar) 

0.01 

kd 

*H,ref 

mol/mpt 

0.005 

*H 

mol/mpt 

0.02 

F 

mol/(mpt  bar) 

9.5 

F 

*CO 

mol/(rapt  bar) 

0.005 

kd 

*CO,ref 

mol/mpt 

5  x  10“15 

F 

K  CO, ref 

mol/mpt 

6  x  10-10 

*0)2,  ref 

mol/(mpt  bar15) 

3.8  x  10-10 

ks 

*CO,ref 

mol/(rapt  bar) 

2.9  x  10“4 

V 

0.1 

UR 

0.5 

aco 

0.15 

Tref 

K 

358 

AH« 

kJmol-1 

9.5 

^  H)o 

kJmol-1 

150 

Eho 

kJ  mol-1 

65 

£co2 

kJ  mol-1 

85 

Eco 

kJ  mol-1 

39.7 

zn 

1 

zco 

2 

The  rate  constants  for  the  electrochemical  oxidation  steps 
are  scaled  with  respect  to  temperature  based  on  the  activation 
energy. 

=  (16) 

We  have  inferred  a  value  of  65  kJmol-1  for  E^0  from  the 
literature  data  on  effect  of  temperature  on  degradation  of  cell 
potential  due  to  CO.  We  currently  do  not  have  a  recommendation 
for  E\ j. 

The  rate  constant  for  the  RWGS  reaction  has  been  scaled 
with  respect  to  temperature  using  85  kJ  mol-1  activation  energy 
(£^q2).  At  present,  we  do  not  have  a  recommendation  for  the 
temperature  dependence  of  k ^  or  k^0.  More  data  on  tempera¬ 
ture  dependence  of  recovery  from  CO  poisoning  by  air  bleed  is 
needed  to  refine  the  choice  of  rate  constants  for  the  oxidation 
reactions  (R-6)  and  (R-7). 

Finally,  catalyst  site  density  (mol-Pt/m2)  is  related  to  the  Pt 
loading  (Lpt)  and  the  specific  electrochemical  surface  area  of  Pt 
(-Apt)  as 

Nc  =  2.2  x  10_5LPtA^,  (17) 

which  is  based  on  2 10  (jlC  cm-2-Pt  as  the  charge  that  Pt  can  hold. 
Gasteiger  et  al.  [10]  recommend  a  value  of  52  m2  g-1  for  Apt 
for  optimized  electrode  structures. 

Figs.  2  and  3  compare  the  modeled  dependence  of  anode 
overpotential  on  CO  concentration  and  temperature  with  the 
single-cell  experimental  data  of  Lee  et  al.  [1]  The  ME  A  con¬ 
sisted  of  Nation  115  membrane  with  0.4  mg  cm-2  Pt  loading 
on  anode  and  cathode  (0.8  mg  cm-2  total  Pt  loading).  The  cell 
was  operated  with  H2/O2  and  H2-CO/O2,  maintaining  1-bar 
partial  pressure  of  H2  and  O2  in  humidified  anode  and  cathode 
streams.  Within  the  range  of  CO  concentrations  analyzed,  Fig.  2 
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0.0  0.5  1.0  1.5 

Current  Density,  A  cm2 

Fig.  2.  Model  validation:  effect  of  CO  on  anode  overpotential  at  85  °C.  Symbols 
denote  experimental  data  from  Ref.  [1]. 

shows  good  agreement  between  the  anode  overpotential  calcu¬ 
lated  using  the  bridge-site  assumption  and  the  experimental  data 
at  85  °C.  Fig.  3  shows  that  for  20-ppm  CO  in  fuel  H2,  the  assump¬ 
tion  of  150 kJ mol-1  for  AHq0  gives  good  agreement  between 
the  calculated  and  measured  anode  overpotentials  over  the  tem¬ 
perature  range  70-100  °C  although  some  disparity  is  observed 
at  40  °C. 

Experimental  data  of  Wilson  et  al.  [8]  was  selected  to  validate 
the  choice  of  rate  constants  for  CO  and  H2  oxidation  reactions 
over  Pt.  These  experiments  were  conducted  with  Nation  105 
membrane  and  a  low  0.14  mg  cm-2  Pt  loading  on  anode  and 
cathode.  The  cell  was  operated  with  humidified  anode  and  cath¬ 
ode  streams  at  80  °C  and  lbar.  Fig.  4  presents  the  modeled 
polarization  curves  for  neat  H2,  5-ppm  CO  in  H2  and  20-ppm 
CO  in  H2,  and  compares  them  with  the  experimental  data.  It 
also  includes  modeled  curves  with  2%  air  bleed  in  H2  contain¬ 
ing  5  ppm  CO  and  curves  with  1  %  O2  bleed  in  H2  containing 
20  ppm  CO.  In  both  cases,  the  choice  of  rate  constants,  k ^  and 
&CQ,  listed  in  Table  1,  is  seen  to  yield  an  excellent  match  with 
the  measured  voltage  recovery  with  air/oxygen  bleed. 

Experimental  data  of  Gu  et  al.  [9]  were  selected  to  validate 
the  CO2  degradation  model.  These  experiments  were  conducted 
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Fig.  3.  Model  validation:  effect  of  temperature  on  anode  overpotential,  20-ppm 
CO  in  fuel  H2.  Symbols  denote  experimental  data  from  Ref.  [1]. 


Fig.  4.  Model  validation:  effect  of  air  bleed  on  anode  overpotential  at  80  °C. 
Symbols  denote  experimental  data  from  Ref.  [8]  (empty  diamonds,  neat  H2; 
empty  squares,  5  ppm  CO  with  2%  air;  empty  triangles,  20  ppm  CO  with  1% 
O2;  filled  squares,  5  ppm  CO;  filled  triangles,  20  ppm  CO). 

using  the  Primea  MEA,  Series  5510  that  has  a  25-pim  thick 
gore-select  membrane  and  catalyst  loadings  of  0.4  mg  cm-2  Pt 
on  both  anode  and  cathode.  The  cell  was  operated  at  70  °C  and 
1  bar,  with  the  anode  and  cathode  gases  humidified  to  80  and 
70  °C  dew  point  temperatures,  respectively.  The  stoichiometry 
was  1.2  for  the  anode  (83%  H2  utilization)  and  2.0  for  cathode 
(50%  O2  utilization).  Fig.  5  compares  the  modeled  and  mea¬ 
sured  anode  overpotential,  defined  as  the  difference  between  the 
cell  voltage  with  neat  H2  and  the  cell  voltage  with  H2-CO2 
mixture  at  the  same  current  density.  Defined  in  this  manner, 
the  anode  overpotential  includes  the  effect  of  CO2  dilution  on 
the  Nernst  potential,  activation  overpotential  due  to  catalyst 
poisoning  by  CO  produced  by  the  RWGS  reaction  and  any  addi¬ 
tional  mass  transfer-related  overpotentials  due  to  the  presence  of 
20-60%  CO2.  The  good  comparison  indicates  that  our  choice  of 
the  rate  constant  for  the  RWGS  reaction  ( krCQ2 )  in  conjunction 
with  the  rate  constants  for  CO  adsorption  and  electrochemi¬ 
cal  oxidation  reactions  can  explain  the  observed  cell  voltage 
degradation. 


Fig.  5.  Model  validation:  effect  of  20-60%  CO2  on  anode  overpotential  at  70  °C. 
Symbols  denote  experimental  data  from  Ref.  [9]. 


R.K.  Ahluwalia,  X.  Wang  /  Journal  of  Power  Sources  180  (2008)  122-131 


127 


4.  Results 

Several  simulations  were  run  using  the  model  described  in 
Section  2  to  parametrically  study  the  effects  of  impurities  on  a 
reference  PEFC  stack  that  generates  90  kWe  gross  (80  kWe  net 
system  power  output)  in  a  pressurized  fuel-cell  system  with  neat 
H2  as  fuel  [2] .  Some  important  design  and  operating  parameters 
and  assumptions  for  the  stacks  are  listed  below. 

(a)  The  cell  voltage  is  0.65  V  at  90  kWe  rated  power. 

(b)  The  cell  is  isothermal;  the  cell  temperature  is  80  °C  at  all 
operating  points. 

(c)  The  pressure  of  fuel  and  air  streams  at  stack  inlet  floats  with 
flow  rate:  2.5  bar  at  rated  flow,  2.0  bar  at  75%  flow,  1.5  bar 
at  50%  flow,  1.3  bar  at  25%  flow  and  1.1  bar  at  10%  flow. 
The  air  management  system  has  a  maximum  turndown  of 
20. 

(d)  The  stack  operates  at  50%  O2  utilization  and  70%  per-pass 
H2  utilization  if  the  air  management  system  can  match  the 
airflow  rate  required  for  a  given  electrical  load.  For  stack 
operation  below  the  maximum  turndown  of  the  air  manage¬ 
ment  system,  the  O2  utilization  is  less  than  50%  and  the  H2 
flow  rate  in  the  anode  channels  is  held  constant.  In  this  oper¬ 
ating  regime,  both  the  O2  and  the  H2  utilization  decrease  as 
the  electrical  load  is  reduced. 

(e)  Both  anode  and  cathode  inlet  gas  streams  are  humidified  to 
60%  relative  humidity  (RH)  at  80  °C  stack  temperature,  i.e., 
the  dew  point  temperature  is  68  °C. 

(f)  The  membrane  electrode  assembly  (ME A)  consists  of  a  25- 
p,m  thick  PFSA  membrane  whose  physical  and  transport 
properties  (density,  water  uptake,  swelling,  ionic  con¬ 
ductivity,  osmotic  drag  and  water  diffusivity)  are  given 
by  correlations  available  for  Nation  112.  The  anode  and 
cathode  catalyst  layers  have  identical  electrode  struc¬ 
tures:  0.4  mg  cm-2  Pt  loading,  Pt/C  =  0.47,  ionomer/C  =  0.8, 
54  m2  g_1-Pt  electrochemical  surface  area,  and  40%  poros¬ 
ity.  It  is  assumed  that  the  kinetics  of  the  electrochemical 
oxygen  reduction  reaction  over  Pt  can  be  adequately 
described  using  the  parameters  (exchange  current  den¬ 
sity,  transfer  coefficients,  etc.)  derived  by  Gasteiger  et  al. 
[10]. 

Fig.  6  shows  the  modeled  polarization  curve  for  a  single  cell 
of  the  stack  with  neat  H2  as  fuel.  It  includes  performance  losses 
due  to  crossover  of  H2,  O2,  and  N2  across  the  membrane.  The 
recycle  ratio  is  125  (0.8%  purge)  to  control  the  buildup  of  N2  in 
the  anode  gas  channels. 

4.1.  Effect  of  recycle 

Fig.  6  also  presents  the  combined  effects  of  CO  impurity  and 
anode  gas  recycle  on  the  polarization  curves.  Two  sets  of  calcu¬ 
lated  results  are  included  in  Fig.  6  (R  =  0  and  100)  to  isolate  the 
effect  of  CO  in  fuel  H2  from  the  effect  of  CO  buildup  in  the  recy¬ 
cled  anode  gas.  The  polarization  curve  for  R  =  0  (once-through 
anode  gas  flow)  indicates  that  1-ppm  CO  in  fuel  H2  decreases  the 
cell  voltage  by  35  mV  at  1.05  A  cm-2.  The  polarization  curve 
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Fig.  6.  Effect  of  anode  gas  recycle  on  voltage  degradation  due  to  CO  in  fuel  H2. 

for  R=  100  (closed  anode  circuit,  1%  purge)  shows  a  further 
decrease  of  25  mV  at  1.05  A  cm-2  because  of  CO  buildup  due 
to  anode  gas  recycle. 

Our  simulations  suggest  that  above  a  certain  threshold  that  is 
determined  by  O2  crossover,  even  a  small  concentration  of  CO 
in  fuel  H2  can  cause  a  measurable  decrease  in  cell  voltage  at  a 
high  enough  recycle. 

Analogous  to  Fig.  6  for  CO,  Fig.  7  presents  the  effects  of 
CO2  impurity  and  anode  gas  recycle  on  the  polarization  curves. 
The  polarization  curve  for  R=  10  shows  that  1%  CO2  in  fuel 
H2  decreases  the  cell  voltage  by  15  mV  at  1.05  A  cm-2.  The 
cell  voltage  decreases  by  50  mV  if  the  recycle  is  increased  to 
100.  Our  simulations  suggest  that  at  small  recycle  ratios,  the 
decrease  in  cell  voltage  is  primarily  due  to  production  of  CO  by 
the  reverse  water-gas  shift  reaction.  At  large  recycle  ratios,  the 
decrease  in  Nernst  potential  due  to  the  lowering  of  the  H2  partial 
pressure  by  buildup  of  CO2  also  comes  into  play. 

4.2.  Degradation  mechanism 

Fig.  8  has  been  constructed  to  clarify  the  modeled  mecha¬ 
nism  of  performance  degradation  due  to  CO  and  CO2  impurities. 
Fig.  8a  displays  the  isotherm  for  H2  adsorption  on  the  Pt  elec- 
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Fig.  7.  Effect  of  anode  gas  recycle  on  voltage  degradation  due  to  CO2  in  fuel 
H2. 


128 


R.K.  Ahluwalia,  X.  Wang  /  Journal  of  Power  Sources  180  (2008)  122-131 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

1.0 

0.8 

0.6 

0.4 

0.2 


(a)  Neat  Hydrogen 


0H 


0M 


(d)  1%C02,  R=10 


(e)  1%C02)  R=100 


Oco 


2ZL 


0.4  0.6  0.8 

Current  Density,  A  cm2 


Fig.  8.  Average  site  coverage  under  different  operating  conditions. 


trocatalyst.  At  low  current  densities,  the  calculated  average  H2 
coverage  ((0i)=J0i  d AM)  is  greater  than  45%.  In  Fig.  8a,  (0u) 
increases  with  current  density  because  of  the  coincident  increase 
in  the  operating  pressure  that  causes  the  adsorption  equilibrium 
in  reaction  (R-l)  to  shift  to  the  right.  At  constant  pressure,  ($h) 
should  decrease  with  increase  in  current  density  because  of  the 
faster  rate  of  the  H2  oxidation  reaction. 

Fig.  8b  (R  =  0)  indicates  that  CO  adsorbs  strongly  on  Pt 
(<@c 0^0-7  at  0.5  A  cm-2  and  higher  current  densities)  and 
blocks  the  access  of  H2  to  the  catalytic  sites.  It  shows  that 
1-ppm  CO  at  80  °C  can  reduce  (On)  to  <0-2  from  >0.6 
with  neat  H2.  Fig.  8c  indicates  that  (0co)  may  approach  0.8 
with  anode  gas  recycle  ( R=  100),  causing  (On)  to  decrease  to 
<0.1. 

Our  model  indicates  that  whereas  HOR  is  confined  near  the 
membrane-catalyst  interface  with  neat  H2,  it  spreads  out  towards 
the  catalyst-gas  diffusion  medium  interface  in  the  presence  of 
CO.  Our  simulations  show  that  the  HOR  becomes  more  uniform 
and  the  catalyst  utilization  improves  as  (0qo)  increases. 

Fig.  8d  and  e  present  the  adsorption  isotherms  in  the  presence 
of  CO2  in  fuel  H2.  CO2  adsorbs  weakly  on  Pt  but  can  produce 
CO  by  the  catalytic  reverse  water-gas  shift  reaction.  For  this 


reason,  the  adsorption  isotherms  with  CO2  impurity  are  qual¬ 
itatively  similar  to  those  with  CO  impurity  except  that  much 
larger  concentrations  of  CO2  and  high  recycle  ratios  are  needed 
to  produce  similar  changes  in  (0qo)  and  (On)- 

4.3.  Optimum  recycle 

A  parametric  study  was  conducted  to  characterize  the  degra¬ 
dation  in  stack  performance  due  to  CO  impurity  in  fuel  H2.  It 
was  found  that  the  stack  efficiency  (ij)  is  a  useful  metric  for 
characterizing  the  performance  degradation.  Here  r)  is  defined 
as  the  ratio  of  the  dc  power  generated  by  the  stack  to  the  lower 
heating  value  (LHV)  of  H2  utilized  in  the  electrochemical  reac¬ 
tion,  H2  that  reacts  in  the  cathode  due  to  crossover  from  the 
anode,  H2  that  reacts  in  the  anode  with  O2  that  crosses  over 
from  the  cathode,  and  H2  that  is  purged  from  the  anode  recy¬ 
cle  loop.  The  results  of  the  study  are  summarized  in  Fig.  9, 
which  shows  that  for  a  given  set  of  operating  conditions,  there 
exists  an  optimum  recycle  (R*)  at  which  the  stack  efficiency  (??*) 
is  maximum.  The  stack  efficiency  is  lower  than  rj*  for  R<R* 
because  of  the  excessive  loss  of  H2  in  purge  gas.  The  stack 
efficiency  is  also  lower  than  for  R<R*  because  of  the  exces¬ 
sive  buildup  of  CO  in  the  anode  channels  and  the  associated 
decrease  in  cell  voltage.  For  the  reference  conditions  of  Section 
3,  R*  is  135  (0.6%  purge)  for  neat  H2,  in  which  case  the  opti¬ 
mum  purge  is  determined  by  the  crossover  of  N2  to  the  anode 
from  the  cathode.  For  impure  fuel  H2,  R*  is  a  function  of  the 
CO  concentration:  the  higher  the  CO  concentration  in  fuel  H2 
the  lower  the  R*.  The  optimum  recycle  ratio  decreases  from  135 
for  neat  H2  to  100  (1%  purge)  if  the  fuel  H2  contains  lppm 
CO. 

Analogous  to  Fig.  9  for  CO,  Fig.  10  presents  the  optimum 
recycle  ratios  as  a  function  of  CO2  impurity  in  fuel  H2.  Unlike 
the  results  in  Fig.  9  which  show  a  gradual  drop-off  in  r|  for  R  >  R* 
and  <  1-ppm  CO  in  fuel  H2,  Fig.  10  shows  a  sharp  decline  in  r\ 
for  R>R*  and  >0.1%  CO2  in  fuel  H2.  As  explained  later,  this 
is  because  large  anode  gas  recycles  at  high  enough  CO2  content 
can  also  lead  to  significant  CO2  buildup  in  the  anode  channel  and 
cause  the  H2  partial  pressure  and,  therefore,  the  Nernst  potential 
to  decrease. 
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Fig.  9.  Optimum  recycle  ratio  as  a  function  of  CO  concentration  at  1 .05  A  cm-2 . 
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Fig.  10.  Optimum  recycle  ratio  as  a  function  of  CO2  concentration  at 
1.05  A  cm-2. 

4.4.  Effect  of  Pt  loading 

Much  effort  is  being  devoted  to  discover  new  alloys  and  elec¬ 
trode  structures  that  allow  lower  Pt  loading  and  reduction  in 
stack  cost.  Although  the  effort  is  primarily  focused  on  the  oxy¬ 
gen  reduction  reaction  (ORR),  reduction  of  Pt  loading  for  the 
hydrogen  oxidation  reaction  (HOR)  is  also  of  interest.  To  this 
end,  a  parametric  study  was  conducted  to  investigate  the  degree 
to  which  reduction  in  Pt  loading  may  adversely  affect  the  toler¬ 
ance  of  HOR  catalyst  to  impurities  in  fuel  H2.  Fig.  11  presents 
some  representative  results  from  a  parametric  study  on  the  effect 
of  CO  impurity  on  efficiency  of  a  stack  with  0.05  mg  cm-2  Pt 
loading.  These  may  be  compared  with  the  results  shown  in  Fig.  9 
for  0.4  mg  cm-2  Pt  in  anode  catalyst  to  ascertain  the  effect  of  Pt 
loading  on  performance  degradation  due  to  CO  impurity.  The 
performance  degradation  due  to  reduced  Pt  loading  is  seen  to  be 
small  if  the  fuel  H2  is  free  of  impurities — a  well  known  result  that 
the  H2  oxidation  reaction  on  Pt  is  facile  and  can  be  conducted  at 
overpotentials  of  <5  mV,  even  at  a  Pt  loading  of  0.05  mg  cm-2. 

Our  simulation  results  suggest  that  under  the  reference  oper¬ 
ating  conditions,  the  adverse  effect  of  CO  in  fuel  H2  is  amplified 
at  low  Pt  loadings.  At  optimum  recycle  ratios,  reducing  the  Pt 


Table  2 


Effect  of  Pt  loading  on  degradation  of  cell  potential  and  stack  efficiency 


Pt  loading 

0.4  mg  cm  2 

0.05  mg  cm  2 

AV 

A 17 

AV 

A  ri 

CO  (ppm) 

0.1 

18.4 

1.5 

42.3 

3.4 

0.2 

28.1 

2.2 

63.3 

5.1 

0.5 

44.1 

3.6 

95.9 

7.8 

1.0 

60.5 

4.9 

124.3 

10.0 

C02  (%) 

0.1 

5.7 

0.7 

4.2 

0.4 

0.2 

7.8 

1.1 

5.6 

0.8 

0.5 

14.3 

1.8 

11.1 

1.4 

1.0 

20.5 

2.5 

17.4 

2.1 

content  of  the  anode  catalyst  from  0.4  to  0.05  mg  cm-2  causes 
an  additional  voltage  drop  of  24-64  mV  if  the  fuel  H2  contains 
0.1-1  ppm  CO  (see  Table  2);  the  corresponding  additional  degra¬ 
dation  in  stack  efficiency  is  1. 9-5.1%.  Also,  the  lower  the  Pt 
content,  the  lower  is  the  optimum  recycle  ratio,  i.e.,  greater  is 
the  purge  loss  of  H2. 

Comparison  of  the  results  in  Figs.  10  and  12  suggests  that 
reduction  of  Pt  loading  in  the  anode  catalyst  may  have  only  a 
small  effect  on  the  degradation  in  stack  efficiency  due  to  CO2 
alone  in  fuel  H2.  A  counterintuitive  result  from  our  simulations 
is  that,  under  the  reference  conditions,  reducing  the  Pt  content 
actually  slows  the  RWGS  reaction  that  produces  the  CO  poison 
so  that  the  cell  voltage  improves  by  1.5-3. 1  mV  (see  Table  2) 
and  the  stack  efficiency  by  0. 3-0.4%. 

4.5.  Buildup  of  impurities 

Fig.  13  presents  the  buildup  of  CO  and  CO2  in  the  anode  gas 
channels  due  to  CO  impurity  in  fuel  H2  when  the  spent  anode 
gas  is  recycled  to  the  anode  inlet  using  the  optimum  recycle 
ratio.  Fig.  13a  shows  the  CO  concentration  at  the  stack  inlet  and 
at  the  stack  outlet  as  a  function  of  CO  concentration  in  fuel  H2 
(all  concentrations  reported  on  dry  basis).  A  number  of  factors 
combine  to  determine  the  CO  concentration  profile  along  the 
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Fig.  11.  Effect  of  CO  on  stack  efficiency  at  reduced  Pt  loading  in  the  anode 
catalyst. 
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Fig.  12.  Effect  of  CO2  on  stack  efficiency  at  reduced  Pt  loading  in  the  anode 
catalyst. 
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Fig.  13.  CO  and  CO2  buildup  in  anode  channels  due  to  CO  in  fuel  H2  at  optimum 
recycle  ratios  shown  in  Fig.  11. 
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Fig.  14.  Enrichment  of  anode  gas  with  CO2  at  optimum  recycle  ratios  shown  in 
Fig.  12  (fuel  H2  contains  only  CO2  impurity). 


anode  channel:  H2  consumption  by  the  electrochemical  reac¬ 
tion,  electrochemical  oxidation  of  CO,  CO  reaction  with  O2  that 
crosses  over  from  the  cathode  channel,  CO  production  by  RWGS 
reaction,  etc.  The  results  in  Fig.  13a  indicate  that  recirculation 
enriches  the  anode  gas  for  reference  stack  operating  conditions 
and  optimum  purge  level.  This  implies  that  CO  conversion  to 
CO2  is  less  than  70%  (H2  utilization)  even  if  fuel  H2  contains 
0.1  ppm  CO  and  the  amount  of  O2  that  crosses  over  is  stoi- 
chiometrically  sufficient  to  oxidize  it  to  CO2.  This  is  because 
of  low  O2  selectivity  for  CO  oxidation  if  the  CO  coverage  is 
small  as  is  the  case  with  0.1 -ppm  CO  concentration.  O2  selec¬ 
tivity  for  CO  oxidation  is  higher  at  higher  CO  concentrations 
but  the  amount  of  O2  crossover  is  stoichiometrically  insufficient 
to  oxidize  it  to  CO2.  In  our  model,  the  primary  mechanism  of 
CO  oxidation  to  CO2  with  1-ppm  CO  in  fuel  H2  is  through  the 
electrochemical  step  (reaction  (R-4)).  Fig.  13a  shows  that  the 
extent  of  enrichment  increases  with  CO  concentration  in  fuel 
H2.  Over  the  0.1-1  ppm  range  of  CO  concentration  in  fuel  H2, 
there  is  45-175%  enrichment  in  CO  at  stack  inlet  and  225-600% 
at  stack  outlet. 

Fig.  13b  shows  that  recycle  can  enrich  the  anode  gas  in  CO2 
to  a  level  that  is  more  than  two  orders-of-magnitude  higher  than 
the  concentration  of  CO  from  which  it  is  produced  by  oxidation. 
We  calculate  that  the  fraction  of  CO  that  is  converted  to  CO2 
is  <50%  if  the  fuel  H2  contains  0. 1-ppm  CO  and  ~20%  if  it 
contains  1-ppm  CO.  At  these  levels  of  CO  conversion,  CO  will 
not  buildup  in  the  recycled  anode  gas  if  the  stack  is  operated  at 
less  than  20-50%  H2  utilization. 

Analogous  to  the  results  in  Fig.  1 3  for  CO,  Fig.  14  presents  the 
buildup  of  CO2  in  the  anode  gas  channels  due  to  CO2  impurity 
in  fuel  H2.  As  in  Fig.  13,  these  results  are  for  reference  stack 
conditions  and  optimum  recycle  ratios  but  with  no  CO  in  fuel 
H2.  Note  that  anode  gas  recycle  can  enrich  the  gas  in  CO2  by  a 
factor  of  10  at  stack  inlet  and  30  at  stack  outlet  (Fig.  14). 

4.6.  Impurity  limits 

A  parametric  study  was  conducted  to  determine  the  decrease 
in  cell  voltage  and  degradation  in  stack  efficiency  as  a  function 
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Fig.  15.  Relationship  between  CO  impurity  in  fuel  H2  and  degradation  in  cell 
voltage  and  stack  efficiency  at  reference  conditions. 

of  impurity  concentration  in  fuel  H2.  Figs.  15  and  16  present 
the  results  of  the  study  for  reference  stack  conditions,  cur¬ 
rent  density  (1.05  A  cm-2)  that  leads  to  0.65-V  cell  voltage 
with  neat  H2  after  accounting  for  N2  crossover,  and  optimum 
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Fig.  16.  Relationship  between  CO2  impurity  in  fuel  H2  and  degradation  in  cell 
voltage  and  stack  efficiency  at  reference  conditions. 
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recycle  ratios.  These  results,  which  may  be  useful  in  devel¬ 
oping  H2  quality  standards,  can  be  interpreted  as  follows.  Say 
that  it  is  desired  to  establish  the  limit  for  CO  impurity  in  fuel 
H2  on  the  basis  that  the  maximum  decrease  in  cell  voltage  at 
1.05  A  cm-2  should  be  restricted  to  10  mV.  Fig.  15  suggests 
that  the  CO  concentration  then  has  to  be  limited  to  <0.06  ppm 
and  that  the  corresponding  degradation  in  cell  efficiency  would 
be  ~0.8%-point.  The  optimum  recycle  ratio  at  this  operating 
condition  would  be  about  100,  see  Figs.  9  and  13  indicates 
that  with  0.06  ppm  CO  in  fuel  H2,  anode  gas  recycle  does  not 
cause  significant  enrichment  in  CO  because  of  its  oxidation 
by  O2  crossover  through  the  25-[xm  membrane.  As  a  point  of 
reference,  Fig.  15  suggests  that  1-ppm  CO  in  fuel  H2  would 
cause  >60  mV  decrease  in  cell  voltage  and  ~5  percentage-point 
degradation  in  stack  efficiency.  Fig.  13  indicates  that  anode  gas 
recycle  with  1-ppm  CO  in  fuel  H2  (optimum  R  <  100  in  Fig.  9) 
builds  up  CO  level  to  >3  ppm  at  stack  inlet  and  ppm  at  stack 
outlet. 

Applying  the  above  criterion  to  CO2,  Fig.  16  suggests  that  the 
CO2  concentration  in  fuel  H2  needs  to  be  <3000  ppm  in  order  to 
limit  the  decrease  in  cell  voltage  to  10  mV  at  1.05  A  cm-2.  The 
corresponding  degradation  in  stack  efficiency  is  ~1.4%-point 
and  the  optimum  recycle  ratio  from  Fig.  10  is  about  30.  Fig.  14 
indicates  that  a  recycle  ratio  of  30  would  enrich  the  anode  gas 
from  2500  ppm  CO2  in  fuel  H2  to  a  CO2  concentration  of  5%  at 
stack  inlet  and  14%  at  stack  outlet. 

Note  that  the  results  in  Figs.  15  and  16  are  applicable  to 
one  set  of  design  and  operating  conditions.  The  decrease  in  cell 
voltage  and  degradation  in  stack  efficiency  due  to  CO  and  CO2 
impurities  depends  on  many  parameters  such  as  the  operating 
pressure  and  temperature,  membrane  thickness  and  Pt  loading 
in  the  anode  electrocatalyst.  The  effects  of  these  parameters  will 
be  explored  in  a  future  publication. 

5.  Conclusions 

We  have  conducted  a  comprehensive  study  on  the  effect  of 
CO  and  CO2  impurities  in  fuel  H2  on  the  performance  of  a 
pressurized  PEFC  stack  for  automotive  applications.  The  major 
conclusions  from  the  analysis  in  response  to  the  questions  posed 
in  Section  1  are  briefly  summarized  below: 

•  Any  CO  impurity  in  the  fuel  adsorbs  strongly  on  the  Pt  catalyst 
and  degrades  the  cell  voltage  by  blocking  the  access  of  H2  to 
the  catalyst  sites.  At  low  CO  concentrations,  <100  ppm,  the 
observed  degradation  can  be  explained  on  the  basis  of  bridge- 
site  adsorption  of  CO,  followed  by  electrochemical  oxidation 
of  CO  to  CO2. 

•  Any  CO2  impurity  in  the  fuel  adsorbs  weakly  on  the  Pt  cat¬ 
alyst.  At  low  CO2  concentrations,  the  observed  degradation 
in  cell  voltage  can  be  explained  on  the  basis  of  CO  produced 
by  the  reverse  water-gas  shift  reaction  between  CO2  and  the 
adsorbed  H2.  At  higher  CO2  concentrations,  the  decrease  in 
cell  potential  due  to  dilution  of  the  H2  also  contributes  to 
performance  degradation. 


•  At  high  H2  utilization,  depending  on  O2  crossover  and  current 
density,  anode  gas  recycle  can  enrich  the  recirculating  gas 
with  CO  impurity.  Recycle  always  leads  to  buildup  of  CO2  in 
the  anode  channels. 

•  The  buildup  of  CO  and  CO2  impurities  can  be  controlled  by 
purging  a  fraction  of  the  spent  anode  gas.  There  is  an  optimum 
purge  level  at  which  the  degradation  in  the  stack  efficiency 
is  the  smallest.  The  stack  efficiency  is  lower  for  purge  higher 
than  the  optimum  because  of  excessive  loss  of  H2  in  purge 
gas.  The  stack  efficiency  is  also  lower  for  purge  less  than  the 
optimum  because  of  the  decrease  in  cell  voltage  due  to  the 
excess  buildup  of  CO  and  CO2. 

•  The  method  and  results  presented  in  this  paper  can  be  used 
to  determine  the  limits  of  CO  and  CO2  impurities  in  fuel  H2. 
The  impurity  limits  are  functions  of  stack  design  and  oper¬ 
ating  conditions  and  allowable  degradation  in  cell  voltage 
and  stack  efficiency.  As  an  example,  Fig.  15  indicates  that 
CO  concentration  in  fuel  H2  should  be  less  than  0.06  ppm 
if  the  decrease  in  cell  voltage  of  the  reference  stack  is  to  be 
restricted  to  10  mV  at  rated  power.  At  this  concentration,  the 
optimum  recycle  ratio  is  about  100  and  the  stack  efficiency 
degrades  by  ~0.8  percentage-point.  Under  the  same  restric¬ 
tion,  the  CO2  concentration  should  be  less  than  2500  ppm  for 
a  corresponding  decrease  of  ~  1.3  percentage-point  decrease 
in  stack  efficiency  at  an  optimum  recycle  ratio  of  ~40. 
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